Possible anomalous Doppler shift effect in superconductor Sr 2 Ru04 
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The effect of the Doppler shift is studied in a model for the a-/3 bands of Sr2Ru04 consisting of 
two hybridized ID bands. Assuming a superconducting gap with nodes in the diagonal directions, 
we examine the oscillation of the surface density of states and the thermal conductivity under a 
rotating magnetic field. Upon varying the strength of the hybridization, the oscillation in these 
quantities is found to exhibit 2D to ID crossover. In the crossover regime, which corresponds to 
the actual S^RuCU, the thermal conductivity exhibits a two-fold-symmetry oscillation, while the 
four-fold-symmetry component in the oscillation is barely detectable. 



A ruthenate superconductor Sr 2 Ru04 1] has attracted 
much attention as a possible candidate for spin-triplet 
superconductor. After a theoretical prediction that the 
pairing should occur in the spin-triplet channel, sev- 
eral experiments 0, 0, have in fact supported this 
possibility. Early predictions were that the mostprob- 
able d- vector is of the form d(k) = (k x ± i/cj,)z.|2(. In 
this case, the gap does not have nodes, so the material 
should be asrapful superconductor. However, recent ex- 
periments [flTE H> H EiJ indicate the existence of nodes 
(or node-like structures) in the energy gap. Stimulated 
by these experiments, several theoretical models with 
gaps having line nodes or node-like structures have been 
proposed El El El El El El On the other hand, ther- 
mal conductivity measurements under rotating magnetic 
field[llll have indicated that the gap is almost isotropic 
within the planes, suggesting indirectly the existence of 
horizontal nodes in the gap, thereby excluding the possi- 
bility of vertical nodes or node- like structures proposed in 
some theoretical studies, [H Eil Eil Ei| These measure- 
ments under magnetic field have motivated our study. 

In the presence of a magnetic field, it is known that 
the energy spectrum of the quasipa rticle is influenced by 
the Doppler shift-ElElEllIllIIl The Doppler shifted 
states around the nodes of the superconducting gap con- 
tributes to the density of states at the Fermi energy. If 
we assume a free-electron-like Fermi surface, the density 
of states is minimized (maximized) when the magnetic 
field is applied parallel to the nodal (antinodal) direc- 
tion, thereby exhibiting a four-fold symmetry oscillation 
upon rotating the magnetic field. This effect can be used 
to probe the direction of the nodes by thermal conduc- 
tivity measurement. [a E ft VM . V2/ A E i El| but since such 
experiments can be strongly affected by phonons, we have 
recently proposed an alternative method : magnctotun- 
neling spectroscopy. [25| Namely, by rotating the mag- 
netic field in the ab plane, the surface density of states 
(SDOS), and thus the tunneling spectra, oscillates, which 
enables us to determine the position of the nodes in the 



gap without using the phase sensitive spectroscopy based 
on the appearance of the Andreev bound states[26|. As 
a case study, we have considered the case of the high 
T c cuprates and an organic superconductor k-(BEDT- 
TTF) 2 X, where we found that the SDOS takes its min- 
imum when the applied magnetic field is parallel to the 
nodal direction, as in systems having free-electron-like 
Fermi surface. It is not at all clear, however, whether 
this tendency holds regardless of the shape of the Fermi 
surface. 

This is exactly where the present study sets in. Here 
we study the effect of the Doppler shift in the a- (3 bands, 
the quasi one-dimensional(lD) bands, of Sr2Ru04. This 
is motivated by some microscopic theories |l4l Il5j which 
proposes the presence of superconducting gaps in these 
bands having nodes or node-like structures in the diag- 
onal direction. Here we concentrate on the a-(3 bands, 
which implicitly assumes that a large nodeless gap opens 
in the 7 band and a small one in the a-[3 bands, so 
that the main contribution to the density of states at 
the Fermi level comes from the latter bands. Consider- 
ing a 2D model in which two ID bands are hybridized, 
and assuming a gap that has nodes in the four diagonal 
directions (which of course corresponds to vertical nodes 
in 3D systems), we calculate the SDOS and the ther- 
mal conductivity upon rotating the direction of the mag- 
netic field. To our surprise, we find that the four-fold- 
symmetry component in the oscillation of these quanti- 
ties, which should be clearly visible in 2D systems hav- 
ing free-electron-like Fermi surface, can be barely seen 
when the hybridization is moderate as in the a- (3 bands 
of Sr 2 Ru04. What is even more striking in this case is 
that the thermal conductivity exhibits a strong two-fold- 
symmetry oscillation reflecting the quasi-one dimensional 
nature of the Fermi surface. 

First let us focus on the tunneling spectrum, namely 
the SDOS. Assuming that the penetration depth is much 
longer than the coherence length, the vector potential can 
be expressed as A(r) = (HXe z ^ x sin#, —H\e z ^ x cos8, 0), 
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FIG. 1: Left: magnetic field H in the 9 direction. Right: the 
2D model considered in the present study. 



where e kxz and e kyz are the dispersion of the d xz and 
d yz bands, respectively (see FigQ. t and // are fixed at 
t = 0.18to and p, = 0.17fo, where to is the unit of the 
energy about leV, while t' is varied as a key parameter 
in the present study which controls the strength of the 
hybridization between the two bands. Appropriate value 
of t' for Sr 2 Ru0 4 should be t' = 0.01 - 0.02i , but t' is 
varied in a wider range to see the crossover between ID 
and 2D. 

As for the gap functions, its absolute value (note that 
only the absolute value enters in eqns.(3)&(4)) are chosen 
as 



where 9 is the angle of the magnetic field measured from 
the x-axis as shown in FigQ and A is the penetration 
length, respectively. In the following, we assume A is 
much larger than the coherence length £ and A(r) can be 
approximated as A(r) ~ A = (HXsind, — HXcos8, 0). 
Then the quasiparticle momenta k x and k y in the x and 
y directions can be given as k x = k x + (H/n£H o )sm0 
and k y — k y — (H/n^H ) cos#, where H = (/> /(7T 2 £A) 
with 0o = h/(2e). Thus, the SDOS at zero energy (the 
Fermi energy) can be expressed as, |2(| 

P%B) = J d^sHsech 2 , (1) 

PsM = \ ^ {l"fc,cr| 2 - E k ^ a ) + 8{UJ - E- k ,cr)] 

k,(T 

+ \vk,*\ 2 [5(u + E k . a ) + 6(lu + £_ felCT )]} . (2) 



E±k,a = ^ 




with rj k ^ = £ fe;(T + £-k,o- and r fejCT = ^jr\ 2 ka +4|A fc ^| 2 , 

with band index er(= a or (3). The two bands a and (3 
result from the hybridization of two ID bands, where the 
energy dispersion £ k ^ a and ^ k ,p are given as follows. 



£k,a = 7;{( e kxz + e k yz) ~ J (^kxz ~ tkyz) 2 + 4g|} (5) 



£fe,/3 = -^Wkxz + Ckyz) + y {tkxz - tkyz) 2 + ^1} 

i-kxz = -2tcos{k x a) - p, £ ky z = -2tcos(k y a) - p, 
g k = —At'sm.(k x a) sin(k y a) 



|A fe , a | = \A k .p\ = A \(cos(k x ) - cos(fc„))|, (6) 

in order to take into account the node-like structures of 
the spin-triplet gap functions (having p x + ip y symmetry) 
found mrefkhjlhj. Although we assume this phcnome- 
logical form for simplicity, our conclusion is not qualita- 
tively affected by the detailed form of the gap, that is, 
only the direction of the nodes is important. In the fol- 
lowing calculation, we fix the parameters Ao = 0.05to 
and T = 0.002io, but our conclusion is not qualitatively 
affected by the choice of these values. 

In FigO we plot the normalized SDOS px = 
p(6,H)/p(0,H = 0.15£/o) as a function of 9 for several 
values of t' . In the case oft' = O.lio, where the hybridiza- 
tion is strong and the Fermi surface is round, px is mini- 
mized at 9 = 7r/4, namely the nodal direction, and maxi- 
mized at 9 = 0, 7r/2, the antinodal direction, which is con- 
sistent with the previous theories for 2D systems. 

H3 By 

contrast, when the hybridization is weak and the bands 
are essentially ID (t' = 0.005£o), the maximum and the 
minimum of the oscilation is entirely reversed. Conse- 
quently, for a moderate hybridization t' = O.Olto, the 
SDOS becomes almost constant upon rotating the mag- 
netic field. 
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FIG. 2: The normalized SDOS plotted as a function of 9 
for H = 0.15.Ho with a: t'/t = 0.1, b: t'/t = 0.01, and c: 
t'/t = 0.005. 

This crossover between large and small t' can be un- 
derstood as follows. Let us first note that the Doppler 
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shift is essentially given by vp • A$ , where Vf is the Fermi 
velocity ^3 an d the vector potential A$ is perpendicu- 
lar to the magnetic field H . When the hybridization is 
sufficiently strong, the Fermi surface is round as in Figs. 
E{al) and[3Ja2), so that the situation is the same as in 
the previous studies for 2D systems. 18] Namely, all the 
states around the diagonal nodes have non-zero vf, x and 
vf, v , thereby contributing to the density of states when 
the magnetic field is applied in the antinodal direction 
(Fig. E{al)). When the field is applied in the nodal di- 
rection, on the other hand, the states around the nodes 
parallel to the field do not contribute to the density of 
states, while the contribution of the states around the 
nodes perpendicular to the field is only a factor of \[2 
larger than it is when the field is applied in the antinodal 
direction (Fig|3fa2)). Thus, the SDOS becomes larger 
when the field is applied in the antinodal direction. 

The situation completely changes when we consider 
the weak hybridization limit, t' — 0, where the two one- 
dimensional bands contribute independently to the den- 
sity of states. In this case, when the magnetic field is 
applied along the antinodal direction (FigEjfbl)&(b3)), 
only the states on the Fermi surface parallel to the field 
gives contribution. On the other hand, for a field in the 
nodal direction, states on both Fermi surfaces have con- 
tribution that is a factor of y/2 smaller than it is when 
the field is in the antinodal direction (FigJS](b2)). Conse- 
quently, in the ID limit, the SDOS becomes larger when 
the field is applied in the nodal direction. 

We now move on to the thermal conductivity. Thermal 
conductivity k xx {8, H) can be expressed as |27ll28j 



,H) 



k,<7 



k,a V k,x T k,<7 



(7) 



with Vk,x = {dEk tCr / dk x ), where we assume a constant 
Tfc.o- with Tfc,<T = tq- As for the energy dispersion E^, we 
assume the form adopted for the calculation of SDOS for 
simplicity. Although this may not accurately correspond 
to the actual experimental situation, we believe that the 
essential physics can be captured within this formalism. 
The thermal conductivity for H = 0.15Hq, normalized 
as kt = k xx (6,H)/k xx (0,H), is shown in Fig. Ufa). For 
t' = O.lto, kt is minimized around 8 = n/A and the 
oscillation essentially has a four-fold symmetry, which is 
consistent with the previous studies for 2D systems. [l3| 
By contrast, for t' = 0.005to and t' = O.Olto, kt ex- 
hibits a strong two-fold-symmetry oscillation, taking its 
maximum at 8 = n/2 and a minimum at 8 = 0. The 
anomalous two-fold-symmetry oscillation in the case of 
weak to moderated hybridization can be understood by 
considering again the t' — limit. Namely, when 8 = 0, 
only the states on the d yz branch of the Fermi surface, 
where Vk x = 0, gives contribution to the density of states 
(see FigOJbl)), so that k xx (8, H) ~ according to eq.0 
at zero temperature. When 8 = n/2, on the other hand, 
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FIG. 3: The states (hatched circles) around the nodes of the 
gap contributing to the density of states when t' is sufficiently 
large and the magnetic field is applied in (al) the antinodal 
or (a2) the nodal directions, or when t' = and the magnetic 
field is applied in the (bl)& (b3) antinodal or (b2) the nodal 
directions. The size of the circles schematically represent the 
magnitude of the contribution. 



the states on the d xz branch contribute to the density of 
states (Fig|3{b3)), so that k xx (8,H) becomes large. 
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FIG. 4: (a)«T and (b)nfT plotted as functions of 6 for H = 
0.15H . a: t'/t = 0.1, b: t' /t = 0.01, and c: t'/t = 0.005. 



Let us now look into the four-fold-symmetry compo- 
nent of k xx (8,H). Since the two-fold-symmetry compo- 
nents in k xx (8, H) and K yy (8, H) have the same absolute 
values and the opposite signs to each other, we can fo- 
cus on the four- fold-symmetry component by looking into 
k fT defined as k/t — Kf{8, H) / Kf (0, H ) 

«/ (^ H ) = \ [«x* (8,H) + K yy (8,H)]. 
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As seen in FigQJb) , the 9 dependence of k /t qualitatively 
resembles that of the SDOS shown in Fig|2] Here again, 
the four-fold-symmetry component in the oscillation of 
k xx (9, H) can be barely seen for moderate (and realistic) 
values of t'. The absence of (or very small) four- fold- 
symmetry component is consistent with the experimental 
observations in refs. but the angles where k xx (6, H) 
is maximized and minimized are reversed compared to 
the two-fold-symmetry oscillation found in ref. |2( . How- 
ever, since we have studied the case of T -C A consid- 
ering that the effect of phonons can be neglected at low 
temperatures in the actual experiments, we believe that 
further experiments at temperatures much lower than the 
energy scale of the superconducting gap is necessary to 
verify our prediction for vertical diagonal nodes. 

To summarize, we have studied the effect of the 
Doppler shift on SDOS and the thermal conductivity in 
a 2D model consisting of hybridized ID bands, assuming 
a superconducting gap that has nodes in the diagonal di- 
rections. When the hybridization is strong and the Fermi 
surface is round, both the SDOS and k xx (8, H) exhibits a 
four-fold-symmetry oscillation, taking a maximum (min- 
imum) when the magnetic field is applied parallel to the 
antinodal (nodal) directions, which is consistent with pre- 
vious theories for 2D systems. By contrast, when the hy- 
bridization is very weak, the SDOS is maximized (min- 
imized) when the magnetic field is applied parallel to 
the nodal (anti- nodal) directions. What is more remark- 
able in the case of weak to moderate hybridization is 
that the thermal conductivity exhibits a strong two-fold- 
symmetry oscillation, reflecting the quasi-lD nature of 
the Fermi surface. Since the actual Sr2Ru04 corresponds 
to the regime with moderate hybridization, k xx (9,H) 
should exhibit a two-fold-symmetry oscillation, taking 
its maximum (minimum) when the magnetic field is ap- 
plied in the y (x) direction, while the four-fold-symmetry 
component in the SDOS and k xx (6,H) should be barely 
detectable. 

Y.T. would like to thank Y. Matsuda for useful and 
fruitful discussions. 
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